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Abstract. We investigate the static and spherically symmetric solutions of Einstein's 
equations for a scalar field with non-canonical kinetic term, assumed to provide both 
the dark matter and dark energy components of the Universe. In particular, we give 
a prescription to obtain solutions (dark halos) whose rotation curve Vc{r) is in good 
agreement with observational data. We show that there exist suitable scalar field 
Lagrangians that allow to describe the cosmological background evolution and the 
static solutions with a single dark fluid. 
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1. Introduction 

The confidence region of cosmological parameters emerging from the analysis of data 
from type la Supernovae (SNIa), Cosmic Microwave Background (CMB) anisotropies 
and the large scale structure of the Universe, suggests that two dark components govern 
the dynamics of the present universe. These components are the Dark Matter (DM), 
responsible for structure formation, and an additional Dark Energy (DE) component 
that drives the cosmic acceleration observed at present. In this paper we focus on Uni- 
fied models of Dark Matter and dark energy (UDM) that can provide an alternative to 
our interpretation of the nature of the dark components of our universe. These models 
have the advantage over the DM + DE models (e.g. ACDM) that one can describe 
the dynamics of the universe with a single dark fluid which triggers the accelerated 
expansion at late times and is also the one which has to cluster in order to produce 
the structures we see today. However, the viability of UDM models strongly depends 
on the value of the effective speed of sound Cs [H El [3] , which has to be small enough 
to allow structure formation [U [5l [6] and to reproduce the observed pattern of CMB 
temperature anisotropies [H El HI [8] . 

Several adiabatic or, equivalently, purely kinetic models have been investigated in 
the literature. For example, the generalized Chaplygin gas ([HI HOl [H] (see also [12j), 
the Scherrer [13] and generalized Scherrer [7] solutions, the single dark perfect fluid with 
a simple 2-parameter barotropic equation of state [l^ , or the homogeneous scalar field 
deduced from the galactic halo space-times [15] . 

Moreover, one can build up scalar field models for which the constraint that the 
Lagrangian is constant along the classical trajectories allows to describe a UDM fluid [7] 
(see also Ref. [I6|, for a different approach). Alternative approaches to the unification 
of DM and DE have been proposed in Ref. [1^ , in the frame of supersymmetry, and in 
Ref. [ISj, in connection with the solution of the strong CP problem. 

One could easily reinterpret UDM models based on a scalar field Lagrangian in 
terms of - generally non-adiabatic - fluids [19], EH] . 

A complete analysis of UDM models should necessarily include the study of 
static solutions of Einstein's field equations. This is complementary to the study of 
cosmological background solutions and would allow to impose further constraints to the 
Lagrangian of UDM models. The authors of Refs. [21] and [15] have studied spherically 
symmetric and static configuration for k-essence models. In particular, they studied 
models where the rotation velocity becomes flat (at least) at large radii of the halo. In 
these models the scalar field pressure is not small compared to the mass-energy density, 
similarly to what found in the study of general fluids in Refs. [22l [23l [2ll [25], and the 
Einstein's equations of motion do not reduce to the equations of Newtonian gravity. 
Further alternative models have been considered, even with a canonical kinetic term in 
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the Lagrangian, that describe dark matter halos in terms of bosonic scalar fields, see 
e.g. Refs. [261 [271 [281 [291 [30]. 

In this paper we assume that our scalar field configurations only depend on the radial 
direction. Three main results are achieved. First, we are able to find a purely kinetic 
Lagrangian which allows simultaneously to provide a fiat rotation curve and to realize 
a unified model of dark matter and dark energy on cosmological scales. Second, we 
have found an invariance property of the expression for the halo rotation curve. This 
allows to find purely kinetic Lagrangians that reproduce the same rotation curves that 
are obtained starting from a given density profile within the standard Cold Dark Matter 
(CDM) paradigm. Finally, we consider a more general class of models with non-purely 
kinetic Lagrangians. In this case we extend to the static and spherically symmetric 
space-time metric the procedure used in Ref. [7j to find UDM solutions in a cosmological 
setting. Such a procedure requires that the Lagrangain is constant along the classical 
trajectories; we are thus able to provide the conditions to obtain reasonable rotation 
curves within a UDM model of the type discussed in Ref. [7]. 

The plan of the paper is as follows. In Section 2 we provide the general framework 
for the study of static spherically symmetric solutions in UDM models. Section 3 is 
devoted to the general analysis of purely kinetic Lagrangians, while in Section 4 we 
analyze more general models with non-canonical kinetic term. Our main conclusions are 
drawn in Section 5. In the Appendix, for completeness we provide the spherical collapse 
top-hat solution for UDM models based on purely kinetic scalar-field Lagrangians, which 
allow to connect the cosmological solutions to the static configurations studied here. 



2. Static solutions in Unified Dark Matter scalar field models 



We consider the following action 



S — Sg + S^p + Sb 



d Xy 



-g 



R 
2" 



£(<^,X) 



+ Sb 



where Si, describes the baryonic matter and 



X 



(1) 



(2) 



We use units such that SnG = c = 1 and signature (— , +, +, +); greek indices run over 
space-time dimensions, while latin indices label spatial coordinates. 
The energy-momentum tensor of the scalar field ip is 



rri<p 



2 6S^ dC{^,X) 



and its equation of motion reads 



dX 



1 



dC 



dC 

dip 



(3) 



(4) 
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We consider a scalar field which is static and spatially inhomogeneous, i.e. such that 
X < 0. In this situation the energy-momentum tensor is not described by a perfect fluid 
and its stress energy-momentum tensor reads 

Tj^u = iP\\ + p)v^ - pg^,u (5) 

where 

P = -P± = , (6) 

= V fj.ip/y/—2X and p\\ = C — 2XdC/dX . In particular, p\\ is the pressure in the 
direction parallel to whereas p± is the pressure in the direction orthogonal to n^. It 
is simpler to work with a new definition of X. Indeed, defining X = —x we have 

n, = W.v/C^xf" (7) 
P\\ = 2x|^ -p. (8) 

Let us consider for simplicity the general static spherically symmetric space-time metric 
i.e. 

ds'^ = - exp (2a (r)) dt^ + exp (2/?(r)) dr^ + r^dn'^ , (9) 

where dfi^ = dO"^ + sin^ 9d(j)'^ and a and j3 are two functions that only depend upon r. 
As the authors of Refs. [211 have shown, it is easy to see that the non-diagonal term 
T^'* vanishes. Therefore (p could be either strictly static or depend only on time. In this 
paper we study the solutions where ip depends on the radius only. 

In the following we will consider some cases where the baryonic content is not 
negligible in the halo. In this case we will assume that most of the baryons are 
concentrated within a radius Vf,. If we define M^, as the entire mass of the baryonic 
component then for r > rf, we can simply assume that M^, is concentrated in the center 
of the halo. 

Considering, therefore, the halo for r > rj,, starting from the Einstein's equations 
and the covariant conservation of the stress-energy (or from the equation of motion of 
the scalar field, Eq. dl])), we obtain 

\{l-[rexp{-2P)]'} = p ^^ = 4npr\ (10) 



dr 

M+M. _|_ py^ 

exp [-2(a + /?)] [r exp (2a)]' - l} = p,, ^ a' = ^^^^^ , (11) 



1 - 



iirr 

[r exp «]'/?'— r(expa)' f=P, (12) 



exp [-{a + 2j3)] 
r 

||--(P„+P) (13) 

(which are the 00, rr and 66 components of Einstein's equations and the r component 
of the continuity equation respectively) where exp (— 2/3(r)) = 1 — (M + M*)/(47rr) and 
R = ln[r^ exp(a(r))]. Here a prime indicates differentiation with respect to the radius 
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r. 



r(expa)' > 0, then we can 
'ore be avoided when trying to 



A first comment is in order here. If i) (3' = and ii) 
immediately see that p < 0. These conditions must theref 
find a reasonable rotation curve. For example, neglecting the baryonic mass, the special 
case of p = A/r^ and exp(a) ~ r™, where A and m are constants, fall into this case. 
We thus recover the no-go theorem derived in Ref. [15] under the assumption that the 
rotation curve <C 1 is constant for all r. 

The value of the circular velocity Vc is determined by the assumption that a massive 
test particle is located at 6' = 7t/2. We define as massive test particle the object that 
sends out a luminous signal to the observer who is considered to be stationary and far 
away from the halo. In this case the value of Vc = rip' can be rewritten as 



2 _ 

1- 



P||rV2 + (M + M,)/(87rr) 



(14) 

P||r72 + 3(M + M,)/(87rr)J 

but when we consider the weak-field limit condition (M + M=K)/(87rr) ^ 1 and since the 
rotation velocities of the halo of a spiral galaxy are typically non-relativistic, fc <^ 1, 
Eq. ([H]) simplifies to [2l] 



A second comment follows from the fact that the pressure is not small compared 
to the mass-energy density. In other words we do not require that general relativity 
reduces to Newtonian gravity (see also Refs. |22l [231 [211 [25]). Notice also that in the 
region where Vc ~ const. <^ 1 it is easy to see that in general exp(a) const, since from 
Eqs. ( ITTj) and (fT5|) one obtains ra' ~ f^. 

Finally, let us point out one of our main results. We can see that the relation (fTSll 
is invariant under the following transformation 

P — ^ P = P + cr(r) p|| — > P|| = P|| + (16) 

if 

3g(r) + rq{r)' = — cr(r) , (17) 

up to a proper choice of some integration constants. Thanks to this transformation we 
can consider an ensemble of solutions that have the same rotation curve. We will come 
back to this point in more detail in the next section. Obviously, these solutions have 
to satisfy the Einstein's equations (fTOi) . (fTTj) and ( |T2l) . and the covariant conservation 
of the stress-energy (fT3!) . Moreover, we will require the validity of the weak energy 
conditions, p > and p\\ + p > 0, i.e. 

r ox 
In the following sections we will consider first a purely kinetic Lagrangian C{X) 
and then two Lagrangians C = f{ip)g{X) and C = g{X) — V{ip). 
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3. Unified Dark Matter models with purely kinetic Lagrangians 

Let us consider a scalar field Lagrangian C with a non-canonical kinetic term that 
depends only on X or x- Moreover, in this section we assume that = (or M ^ M*). 

First of all we must impose that C is negative when X < 0, so that the energy 
density is positive. Therefore, we define a new positive function 

gsix) ^ -C{X) . (19) 

As already explained in Ref . |2lj , when the equation of state p\\ = p\\ (p) is known, one 
can write the purely kinetic Lagrangian that describes this dark fluid with the help of 
Eqs. ([6]) and (IHl). Alternatively, using ( fT3l) . one can connect py and p in terms of r 
through the variable R. Moreover, it is easy to see that starting from the field equation 
of motion (jll), there exists another relation that connects x (i-S- X) with r. This relation 
is 

dgsix) 
L dx 

with k a positive constant. If we add an additive constant to Qsix)^ the solution fl20l) 
remains unchanged. One can see this also through the Eq. ( fT3|l . Indeed, using Eqs. ([6]) 
and ([8]) one immediately finds that Eq. (fT3|) is invariant under the transformation 
p p + K p\\ ^ p\\ — K . In this way we can add the cosmological constant i^' = A to 
the Lagrangian and we can describe the dark matter and the cosmological constant-like 
dark energy as a single dark fluid i.e. as Unified Dark Matter (UDM). 

Let us notice that one can adopt two approaches to find reasonable rotation curves 
fc(r). A static solution can be studied in two possible ways: 

i) The first approach consists simply in adopting directly a Langrangian that provides 
a viable cosmological UDM model and exploring what are the conditions under 
which it can give a static solution with a rotation curve that is fiat at large radii. 
This prescription has been already applied, for example, in Ref. [2T] . 

ii) A second approach consists in exploiting the invariance property of Eq. f|T5l) . with 
respect to the transformation ( |T6i) (when the condition ( ITTj) is satisfied). Usually in 
the literature one reduces the problem to the Newtonian gravity limit, because 
one makes use of a CDM density profile, i.e. one assumes that in Eq. (fTSll . 
p\\ <C M/(47rr^). We can therefore use Eqs. (fT6ll and (fTTll to obtain energy density 
and pressure profiles p{r) and p\\{r) that reproduce the same rotation curve in a 
model with non-negligible pressure. Next, we find an acceptable equation of state 
P\\ = P\\{p) such that we can reconstruct, through Eqs. and ([8]), the expression 
for the Lagrangian C Such a procedure establishes a mapping between UDM and 
CDM solutions that predict the same halo rotation curve Vc{r). As a starting point 
we could, of course, use very different CDM density profiles to this aim, such as 
the modified isothermal-law profile [HI], the Burkert profile [32], the Moore profile 
[33] . the Navarro- Frenk- White profile [Ml [35] or the profile proposed by Salucci et 
al. (see for example [36]). 



k 



exp a(r)l" 



(20) 



Halos of Unified Dark Matter Scalar Field 



7 



As we have already mentioned, the possible solutions one finds in this way have to 
satisfy the Einstein equations ffTOj) . f|TT]) and f|T2l) . the conservation of stress-energy 
(fT3|) and the weak energy conditions. Moreover, the resulting UDM scalar field 
Lagrangian must be able to provide cosmological solutions that yield an acceptable 
description of the cosmological background (see, e.g., Ref. [7]) and low effective 
speed of sound (see for example Refs. [3l |6]) so that cosmic structure formation 
successfully takes place and CMB anisotropies fit the observed pattern ^ [5], |8] . 

Below, using approach i), we provide a worked example of a UDM model with purely 
kinetic Lagrangian which is able to describe a fiat halo rotation curve and then, using 
approach ii), we give a general systematic procedure to obtain a possible Lagrangian of 
UDM model starting from a given CDM density profile. 

3.1. Approach i): The generalized Scherrer solution 

Let us consider the generalized Scherrer solution models obtained in Ref. These 
models are described by the following Lagrangian 

C = -A + gn{X-Xy (21) 

where > is a suitable constant and n > 1. The case n = 2 corresponds to the 
unified model proposed by Scherrer [13]. If we impose that today [{X — X)/X]" <^ 1, 
the background energy density can be written as 

p(a(t)) = Pa + pdm , (22) 

where p\ behaves like a "dark energy" component {p\ = const.) and pdm behaves like 

a "dark matter" component i.e. pdm o,"^, with a{t) the scale factor. 

A static solution for the generalized Scherrer model can be obtained in two possible 

ways: 

1) Starting from the analysis of Ref. [15], in the case of a barotropic Lagrangian for 
the homogeneous field. The authors of Ref. [15j indeed concluded that for n 3> 1 
fiat halo rotation curves can be obtained. In particular they studied spherically 
symmetric solutions with the following metric, 

ds^ = -(^-^^^ dt^ + N{r)dr^ + r^dn^. (23) 

where is a suitable length-scale and / = 2f^. In the trivial case where N{r) 
is constant they find C{X) oc X^^' with / -C 1. For X ^ X the Lagrangian 
C = —A + gn{X — X)^ takes precisely this form. 

2) In the analysis of Ref. [21], solutions where (f is only a function of the radius are 
considered. When the Lagrangian has the form C oc X", with n ~ 10^ the halo 
rotation curve becomes fiat at large radii. In this case n must be an odd natural 
number, such that the energy density is positive. Our model is able to reproduce 
this situation when the matter density is large, i.e. when \X\ ^ X. 



Halos of Unified Dark Matter Scalar Field 



8 



Alternatively, if we wish to avoid large n (c.f. case 2) above) we can start from the 
following Lagrangian 

C = -A + exgn{\X\-Xy (24) 

where ex is some different iable function of X that is 1 when X > X and —1 when 
X < —X < 0. In this way when X > X > we recover the Lagrangian of the 
generalized Scherrer solutions. When X < and x = ~X > X we get 

C=-A-gn{x-Xy (25) 

and, with the help of Eqs. ([6]) and (IE]), we obtain 

P = -P± = -C, = i2n - 1)(7„ {x-Xy + '^ngr.X (x - x)""' - A . (26) 



Now, requiring that x be close to X (i.e. (^x ~ X^ X) and 2ngnX (^x ~ X^ ^ 
0(A), and starting from the relation fl20|) that connects x with r, we get 

^ ~ ng^X^^ r2exp(a(r)) ' ^^^^ 
Consistency with our approximations implies that we have to consider the following 
expressions for radial configurations with r bigger than a minimum radius r^m- In this 
case p\\ and p become 

A B 

P = r„...,^...„^M„/r.-n (28) 



r2exp(a(r)) ' [r^ exp (a(r))]"/("' ^) 

r -\n/(n—l) 

where A = 2{kXf'^ and B = g^ [k^'^ /{ngnX^/'^)\ 

Using the Eqs. ffTUl) and ffTTl) . we be able to calculate the values of the metric terms 

exp (a) and exp {(3) and, thus the value of p and p\\. Alternatively we know that when 

Vc ~ const. <^ 1 at large radii, in a first approximation, we can set exp(a(r)) C = 

const. Therefore for 7^ 3, we can write the function M as 

, AtiB fn — l n-s \ 
M{r) ^ -^-jj—s —r^^ + D (29) 



Cn/{n-l) Vn-3 

where we could also set = for n > 3. Instead, when 1 < n < 3, the second term 
has to be larger than the first one. 
In these cases w ^ becomes 

2/ N ^ , B fn-l 1 , D\ 



For n = 3 we have 

(73/2 



AttB 

M{r)^-— \n [■-_]+ M{f) (31) 



where r > r and 



9 ^ 5 1 , M(f) , , 

-2«^ + ^r'nbU^- (32) 



r vr/ svrr 

In other words we see that the circular velocity becomes approximately constant for 
sufficiently large r. 
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However, let us stress that exp(a(r)) cannot be strictly constant, and that it should 
be chosen in such a way that the positivity of Eq. f|T2|) is ensured. 

This example can be generahzed also to M^, 7^ 0. Obviously, in such a case we have 
to assume that r > rf, > r^m. In this case k, rmin, A, B (through exp(/3(r))) and C 
depend on M*. 

The spherical top-hat solution for this model, which provides the link with the 
cosmological initial conditions, is described in the Appendix. 

3.2. Approach ii): A general prescription to obtain UDM Lagrangians starting from a 
profile of an energy density distribution of CDM 

Defining the energy density distribution of CDM as pcDuir) (with pcdm = 0), the 
transformation f|T6l) becomes 

p{r) = pcDuir) + cr{r) , p^^{r) = q{r) . (33) 

Now, starting from a given CDM density profile, through Eqs. ( |TOl) . ( ITTi) . ( |T3l) and 
( IT7|) we can determine exp (a) , exp (/3), p and p\\ . In a second step we provide the 
conditions to ensure that the energy density is positive 0. In this case, after some simple 
but lenghty calculations, we find 

2V) (r^^^^^-2rQ(r))-2Q^(r) 

, ni \ l A , o ^CDM(r) 3 \ rMcDM(r) ^ 2 \ /g.N 
+ Q(r) I 4r + 3 — h 4r pcDM I = — (4 + 3r pcDM j , (34) 

B[r) = Q{r) , (35) 

1 - Q'ir) 

a{r) = (37) 

where Q{r) = r{r'^q + 1) , B{r) = rexp (—2/3) and A{r) = {ra' + 1). Here we define 
^cdm('^) = 47r Jq r"^ pcDuif^) dr. At this point it is easy to see that Eq. (IMl) does not 
admit a simple analytical solution for a generic Pcdm- On the other hand we know that, 
through pcDM; all these functions depend on the velocity rotation curve Vc{r). Moreover 
v'^{r) <^ 1. Therefore, defining Vc as the value that Vc assumes when the rotation curve 
is flat at large radii or the maximum value of Vc with a particular profile of Pcdm, we 
can expand Q, A and B as 

Q(r) = Q(o)(r) + vlQ^,^{r) + ^Q(2)(r) + . . . , 

Air) = ^(o)(r) +t;,M(i)(r) + Sl!^A(2){r) + . . . , 

I Thanks to this condition, through Einstein's Eq. (|12|) . we can evade the no-go theorem derived in 
Ref. [S]. 
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B{r) = S(o) (r) + v'^B^,^ (r ) + ^^(2) (r) + . . . . (38) 
Following this procedure we can determine p and p\\ in a perturbative way, i.e. 

P(^) = P(o) (r) + (r) + ^p(2) (r) + . . . , (39) 

Pll (^) = Pll (0) (r) + (1) (r ) + ^pii (2) (r ) + . . . . (40) 

Now, looking at the various CDM density profiles which have been proposed in the 
literature [211 E21 ESI EH ESI ES] , we see that we can always take Pcdm as 

PcDM(r) = WcPcDM (i)(r) , (41) 

then 

rr 

McDM(r) = v^^McDM (i)(r) = An f^CDM (i)(r) df . (42) 
For the zeroth-order terms we immediately obtain 
2(0) = r , 
-4(0) = 1 , 

i3(o) = r . (43) 
At the first order we get 

(f) dr , 



2(1) = 




r -3 


r . 


r Pcdm (i) 


Ai) = 


1 


McDM (l)('^) 


2r 




%) = 


2 


rr 


r J 


1 ^VcDM(l) 




. _ ^CDM(i)(r) _ 

• ^ 47r ^ ^ 



For completeness we write also the second order for Q 



^ 1 r McDM(l)(f) [2 ,2 

2(2) = c^r — ^-Q(y(r) -r pcDM(i)(r 



(45) 



Let us stress that if one considers also terms O (v^), Eq. (fn|) instead of Eq. (IT^ should 
be used. In such a case, fc slightly changes with respect to the velocity rotation curve 
that one obtains using a CDM density profile. 

For our purposes we can consider only the zeroth and the first-order terms. At this 
point, we can finally calculate the value of p and . We get 

p{r) = pcDM(r) + ^ = v"^ (^^ rVcDM (i)(r) df - pcoM (i)(r) j , (46) 



Q(r) - r _2 2 



As far as the values of the metric terms exp (a) and exp (/3) are concerned, we obtain 
the following expressions 

"_2 r 1 McDM(l)(r) 



exp (2a) = exp (2a (r)) exp 



r 47r 



(4^ 
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exp (-2/3) = 1 + ^ (2 1 r-pcDM (1) (f ) c^r - r^^^^^^^ ] . (49) 



Now, it is immediate to see that if we want a positive energy density we have to 
impose 2 Jq r^pcDM df > r'^ Pcdm{i)('")- From Eq. (fTOj) we know that M(r) = 
47r //j^j fV(r) +M(f(o)) and Mcnuir) = 4:1t pcBuif) dr + McDuir). Therefore we 
need to know what is the relation between f and f(o). This condition is easily obtained 
if we make use of Eq. f|T5l) . Indeed, we get 

M(i)(f(o)) -McDM(i)(r) ^ 2 rm ^3 ... 

^^^^ — ^ — h — / rycDM{i){r) dr = r pcDM(i)(r) dr , (50) 

47r r(o) JO Jf 

which finally guarantees the invariance of the rotation velocity with respect to the 
transformation in Eqs. ( JT6l) and (JT71) . 

Let us, to a first approximation, parametrize the various CDM density profiles, at 
very large radii (i.e. when we can completely neglect the baryonic component) as 

PCDM = (51) 

where k is a proper positive constant which depends on the particular profile that is 
chosen [31], [32], [SS] [3H |35l |36] . For example for many of the density profiles the slope 
is n = 3 for large radii [321 [331 [31 [35l [36] . 

In this case a positive energy density p > requires n > 2. At this point let us 
focus on the case where 2 < n < 4, since this gives rise to the typical slope of most of 
the density profiles studied in the literature. Therefore we obtain for p(r) and p\\{r): 

pir) = vl — -— , Pll(^) = T— — • (52) 

In particular, 

1) for 77, = 2, we get 

p(r) = 0, p||(r) = pcDM =vl , (53) 

and for the relation between f(o) and f one can choose, for example, f(o) = f = 0. 
In other words, for large radii we have that p(r) <^ P||('")- 

2) Also for 2 < 77, < 3 one can choose f(o) = f = 0. 

3) For n = 3 

P{r) = pcDM, P\\ {r)=vl , (54) 

and, through Eq. f [50l) . we have to impose that 

M(i)(r(o)) -McDM (i)(r) ^ _ ^ ^^^^ 

Notice that the energy density profile is the same as the CDM one only for large 
radii so that M(i)(r) differs from Mcdm (i)('^)- 



Halos of Unified Dark Matter Scalar Field 



12 



4) In addition, for 3 < n < 4, also through Eq. (jSU]), we have to impose that 

M(i)(f(o))-McDM(i)(r) _ f^-" in - 2) 

47r n-3 (4-n)(3-n) ^ ' ^ ' 

Now let us focus where 2 < n < 4. Starting from Eq. ( l52l) to express = p\\{p) we 
solve Eq. ([8]) to recover the Lagrangian for the scalar field 

2k . 

p(x) = -£ = /cx , Pix) = ^n_2) ^^^ ^^^^ 

where is a positive integration constant. We can see that, for this range of n, the 
exponent is larger than 1; thus there are no problems with a possible instability of the 
Lagrangian (see Refs. [211 EH]). Therefore, through the transformation p ^ p + A 
P\\ ~^ P\\ this Lagrangian can be extended to describe a unified model of dark matter 
and dark energy. Indeed, starting from the Lagrangian of the type (12^ . when |X| ^ X 
and if A; = gn, C takes precisely the form fl57|) . 

Finally, we want to stress that this prescription does not apply only to the case 
of an adiabatic fluid, such as the one provided by scalar field with a purely kinetic 
Lagrangian, but it can be also used for more general Lagrangians C{ip,X). 



4. Unified Dark Matter models with non-purely kinetic Lagrangians 

Let us consider more general Lagrangians of type C = C{(f,X), with a non-canonical 
kinetic term, in order to find a UDM model with acceptable cosmological speed of sound. 
In this case we have one more degree of freedom: the scalar field configuration itself. 
Therefore, we have to impose a new condition to the solutions of the equation of motion. 
Ref. [7] required that the Lagrangian of the scalar field is constant along the classical 
trajectories. We want to know whether such a condition could be applied to the static 
spherically symmetric space-time metric. We would also like to know what the behavior 
of the rotation velocity Vc in the halo of a spiral galaxy is like for this class of models. 
In the next subsections we will consider first a Lagrangian of the form C = f{ip)g{X) 
and then a Lagrangian of the form C = g{X) — V{(f). In these cases, for simplicity, we 
will assume that f{ip) and V{ip) are positive. 



4.I. Lagrangian of the type C = f{(f)g{X) 



Let us write the Lagrangian in the form C = f{{p)g{X) = —f{^p)gs{,x)- Immediately we 
notice that the requirement of having a positive energy density imposes that gs{x) is 
positive. In this particular case the equation of motion (jl]) becomes 



_d_ 
dR 



In 



2x^-9s{x) 



+ 2x 



dx 



2X 



dx 



H -1 



gsix) 



dlnfi^) 
dR 



Moreover from Eq. ([8]) we obtain for p\\ 

d^n [gsix)] 



fiv)gsix) 2x- 



dx 



i5t 



(59) 
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Following the procedure previously explained we impose the constraint C = —p = —A, 
i.e. 

fi^) = ^ , (60) 
9s{X) 

which, inserted in the equation of motion fl58|) . allows to find the following general 
solution 

/^^[9six)] ^ k/2 ^ 
dx r'^ exp (a) ' 



where is a constant of integration. Now, inserting Eqs. mOh and (1611) into the relation 
(|59|1 . we obtain 

Ak 

PW = -2 - ^ • (62) 

r'' exp [a) 

Using this expression and considering the halo for r > and M ^ M^,, we are finally 
able to get the expression for the rotation velocity 

2 Ak/2 Ar2 ^ ^ 

vl ^ ^ - ^ • 63 

exp [a) 3 

If exp (a) ~ const, this expression leads to a flat rotation curve for all radii r < rmax 
such that r^jjx ^ 3k/ {2 exp (a)) and provided that the constant k is positive. Therefore, 
in the future we will always neglect the second term in Eq. ( l63l) . 

It is important to stress that the results outlined in Eqs. (I60l)-(l63l) give an efficient 
recipe to obtain a fiat halo rotation curve within the UDM scenario. Once a Langrangian 
(i.e. g{X)) leading to a viable UDM model on cosmological scales is obtained by 
imposing the constraint C = —A (see Ref. [Z]), a fiat rotation curve is guaranteed through 
Eq.( l63i) . There are however two important requirements that have to be satisfied. The 
function Qsix) must allow for a positive integration constant k through Eq. (!6T!) . and the 
Lagrangian must satisfy the stability conditions discussed in Refs. [211 EH EH]), which 
require dC/dX > and dC/dX + 2Xd'^C/dX'^ > (so that the speed of sound is 
positive both in the cosmological setting and for the static solution). 

In the second part of this subsection we will consider first a situation where M* = 0, 
in other words when our halo is composed only of the dark fiuid, and then a situation 
where there is a non-negligible baryon contribution in the inner part of the halo. 

4.I.I. Case = 0: halo composed only of the dark fluid Starting from p = A and 
p\\ = AA;/(r^exp (a)) — A we can exphcitly calculate the value of exp (a) and exp (/3) 
through Eqs. ffTOl) and ffTTl) . Therefore, for M^, = we get 

Ar^ 

exp i-2/3) = 1 - — , (64) 



exp (a) = — 



Ak I / Ar2\^/^ 



1/2 

r 



AA; 1 - Ar^V/^ 



3 



+ 1 , (65) 
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where k is a suitable positive integration constant. In particular, the value of k should 
be such that the term on the RHS of Eq. fl65l) is positive, i.e. 

-1 

(66) 



AX 1/2 
-] r> 



cosh TT + 1 

VAA; 



It is very important to stress that, in this case, the weak energy conditions are satisfied. 
In other words, through this prescription, we are able to evade the no-go theorem derived 
in Ref. [15]. 

Using Eq. ( |65i) and Eq. (1631) we can obtain the following expression for the circular 
velocity 



Voir) 



1/2 



2k , 



(t) 



1/2 



r 



1-1 



3 



+ 1 



(67) 



In order to have values of fc ~ 10 ^ we must impose that 2k/ (Ak) ~ 10^ ^ 3/(Ar^„^.). 
Imposing this condition we can obtain an approximately flat halo rotation curve. 

A simple inspection of Eqs. ( !65l) . (166|) and (1671) shows an interesting property of our 
result. There is a minimum radius ~ (A/3)^^/^ [cosh (2/t/(AA;))]~^ required for the 
validity of fl66l) . Obviously it is necessary that r^m ^ rgai (<^ rmax) where Vgai is the 
typical radius of our halo. 



4-1-2. Case ^ 0.' non-negligible baryonic component in the center of the halo 

In this subsection we assume that r > r^. If M ^ M^, we recover the same result of 
the previous subsection; if ^ O(Ar^), using Eqs. (fTUjl and (ITT]) , we obtain 



exp (-2(3) 



Attt 



exp [a) 



Ak 



1 - 



Arrr J 



— ^ + cosh 
Afc 



47rr\V2 



(6^ 



(69) 



where suitable positive integration constant. In particular it easy to see that 

and k (through exp (/?(r))) depend also on the value of M^,, since one is considering 
r > Vb. Obviously, these functions exist only for r > r* > M*/(47r), having defined r^: as 
the value of the radius for which exp (a(r^,)) = 0. In this case, using the approximate 
relation ( fT5l) . Vc reads 




AirrJ 



+ cosh" 



1/2- 



-1 



1 



(70) 



AA; ' 

To have halo rotation velocities Vc ~ 10"'^ for r ^ r*, we need to impose 2k*/ (Ak) ~ 10^. 
One can see that this condition leads to r^, ^ M*/(47r). Moreover, also in this case we 
have a minimum radius rmin such that fc('^mm) = 1- Starting from Eq. f|T^ we get 
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4-2. Lagrangian of the type C = g{X) — V{ip) 



In this subsection we briefly discuss Lagrangians of the type C = g{X) — V{ip). Let us 
rewrite £ as £ = —[gsix) + ^(v^)]- order to have p > we impose that gsix) > 0. 
For these Lagrangians the equation of motion (jl]) becomes 

d\nV{^) 



X- 



dx 



d In 



X 



dR 



+ 2 



dR 



(72) 



Requiring that the Lagrangian of the scalar field is constant along the classical trajectory, 
i.e 



Vi^) = -gsix)+^, 
from the Eq. ( I72l) we get 

dgsix) k/2 



(73) 



(74) 



dx r"^ exp a 
where /c is a positive constant. 

Now, inserting Eq. (!74l) into Eq. ([8]) we obtain the same expressions for and Vc that 
we obtained in the last subsection i.e. Eqs. (1621) and (l63l) . respectively. 



5. Conclusions 



In this paper we have investigated static spherically symmetric solutions ( "dark halos" ) 
of Einstein's equations for a scalar field with non-canonical kinetic term. Assuming 
that the scalar field depends only on the radius, we studied Unified Dark Matter models 
with purely kinetic Lagrangians. In particular, we obtained a purely kinetic Lagrangian 
which allows simultaneously to produce fiat halo rotation curves and to realize a unified 
model of dark matter and dark energy on cosmological scales. Moreover, we gave a 
prescription to obtain UDM model solutions that have the same rotation curve Vc{r) 
as a CDM model with a specified density profile. Next, we considered a more general 
class of Lagrangians with non-canonical kinetic term. In this case we have one more 
degree of freedom (the scalar field configuration itself) and we need to impose one more 
constraint. To this aim, we required that the Lagrangian is constant, C = —A along 
the solutions of the equation of motion. We have studied whether this condition can be 
applied to the static spherically symmetric space-time metric and what the behavior of 
Vc{r) is for this class of models. 

Let us finally stress that these solutions allow for the possibility to find suitable 
Lagrangians that describe with a single fluid viable cosmological and static solutions. 
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Appendix A. Spherical collapse for purely kinetic scalar-field Lagrangians 

Let us assume a flat, homogeneous Friedmann-Robertson- Walker background metric i.e. 

ds"^ = -dt^ + a{tf6ijdx'dx^ , (A.l) 

where a{t) is the scale factor, 6ij denotes the unit tensor. Moreover, for this particular 
case, the Hubble parameter H is a function only of the UDM fluid = p/3. 
Now let us consider a top-hat spherical overdensity with the purely kinetic model with 
the Lagrangian C = —A + gn{X — X)" and with Qn > [7\. For this particular case 
within the overdense region we have a single dark fluid ondergoing spherical collapse, 
which is described by the following equation 

^ = -l{pR + 3pR) (A.2) 

where R(t), and pn are respectively the scale-factor, pressure and energy density of 
the overdense region and where the dot denotes differentiation w.r.t. the cosmic time t. 
Now, pji and pr are defined by the following expressions [7] 

{2n-l)gn{XR-Xr 



PR = A + 2ngnX{XR - X^-' 

PR = gR = -A + QniXR - X)" 

with Xr = X{R) a function of time. 
The equation of motion is 



(A.3) 
(A.4) 



dg 



R 



dX 



+ 2X 



R 



d'^gR\ (DCr 
dN, 



3 2X 



R 



dg 



R 



R 



dX 



. 



R, 



where dNR = dR/R. The solution of Eq. (lA5ll (for dgR/dXR,XR ^ ) is 

2 



X 



R 



dgi 



\dXRj 

where we can choose kR 



kuR- 



(A.5) 



(A.6) 



Xr(^S' 



\dXRj 

Replacing Eq. (1A.4I) in Eq. (1A.6I) we find 

ngn{XR-Xr~'^' 



with Rfn the value of R at turnaround. 



ta 



X 



R 



kRR-^ (A.7) 

Using now the explicit expressions for pr and pr we arrive at the following set of 
equations 

R__l 
R~~3 

2{n-l) _ ^R 



A + ng^XiXR - X)^-' + {n + - Xy 



(XR-xy^-' + xiXR-x) 

For {Xr - X)/X < 1 Eq. (lASl) becomes 
R 1 



ri^gl 



R- 



(A.8) 
(A.9) 



R 



A + ng^\XR,^-Xr\XR,X)^^ 



(AlO) 
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We can now write all the equations that describe the spherical collapse 

2 1 



a 
a 

PA 

Pdm 

R 
R 

PRom 



2 (pa + Pdm) 
A 



2ngn\Xta-Xr'{XtaX)-2 [ — 



1 / a 



g (.PRdm 



2ng^\Xta-Xr-\Xii,J)^^ 



R 

Rta 



(A.ll) 
(A.12) 
(A.13) 

(A. 14) 
(A.15) 



where = a{tta). 

Following now the same procedure of Ref. [39] we can define x and y 
a 



X = 



y = 



a-ta 
R 

'r' 



ta 



In this way we can redefine pdm and Pr^^ such that 



Pdm 



P/?E 



c 



X'^ 



(A.16) 
(A.17) 

(A.18) 
(A.19) 



where JIdm is the (k-essence) dark matter density parameter, and C = (p/pdm)|x=i- 
Then Eqs. flA.llI) and flA.141) become 



dx 
dr 

fy 

dr^ 



(x^DMix)) 2 
1 



1 



1-Q 



DAI 



'x = V. 



^DAlix 



-X 



(A.20) 

(A.21) 

(A.22) 



where dr = Hta\I^DM{x = 1) and Ka = pa/[3H^^^Idm{x = 1)]. 

Defining U as the potential energy of the overdensity and using energy conservation 
between virialization and turnaround, 



we obtain 



where 



rr 

^+2M 



;i + q)y - 2gr 



ta : 



y=l 



C 



(A.23) 
(A.24) 
(A.25) 



in full agreement with Ref. 
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